Paramaterizations of inclusive cross sections for pion production in proton-proton 

collisions. II. Comparison to new data. 
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A set of new, precise data have recently been made available by the NA49 collaboration for 
charged pion production in proton-proton and proton-Carbon reactions at 158 GeV. The current 
paper compares this new data to five currently available arithmetic parameterizations. Although a 
precise fit is not expected, two of the parameterizations do not work very well but the other three 
are able to provide a moderately good, but not precise fit to the proton-proton data. The best two of 
these three parameterizations are scaled to the proton-Carbon data and again provide a moderately 
good, but not precise fit. 
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I. INTRODUCTION 

The NA49 collaboration H [2] has recently completed 
a series of measurements for charged pion production in 
proton-proton (jpp) and proton-Carbon (pC) collisions 
at a beam momentum of 158 GeV, corresponding to 
^fs =17 GeV. This surpasses previous data in that the 
new data is of much higher precision and quality and 
can therefore be used to provide more precise tests of 
hadronic production models. The NA49 collaboration 
[1 [IJ indicate that simple arithmetic parameterizations 
are unable to account for the fine structure seen in their 
data, and they therefore provide a numerical interpola- 
tion scheme. We agree with this. However arithmetic 
parameterizations are used in a wide variety of applica- 
tions including simulation of particle physics experiments 
Q , simulations of cosmic rays showers in the Earth atmo- 
sphere H [H, 0| , description of particle reactions relevant 
to astrophysics Jjand predicting radiation environments 
inside spacecraft [8|. For those using such parameteriza- 
tions it is of interest to know how they compare to the 
new precision NA49 data [J,[3, even though a precision 
fit will not be possible. Blattnig, Swaminathan, Kruger, 
Ngom, and Norbury Q analysed a set of parameteriza- 
tions currently available and compared to an extensive 
data set for both neutral and charged pion production in 
pp collisions. It was concluded that the parameterization 
of Badhwar and Stephens [l(| provided the best overall 
description of charged pion production. It is interesting 
to see how this compares to the new data [l|. 

The Blattnig et al parameterizations @ have also been 
used in a variety of astrophysical and astroparticle appli- 
cations |I1 [11 IS Q [la Hill where the interest is in 
calculating the spectrum of gamma rays, electrons and 



neutrinos which result from the decay of pions produced 
in proton-proton interactions. For example, Bernado [Til ] 
used the high energy pp —> ir° cross sections in order to 
calculate the 7-ray spectrum from microquasars. This 
spectrum can be measured using gamma ray telescopes 
on satellites. Pion production cross section parameter- 
izations @ have also been used recently in nuclear and 
particle physics applications [H, [H H3, UK- The nu- 
clear modification factor Raa is basically the ratio of a 
nucleus-nucleus cross section divided by a scaled nucleon- 
nucleon cross section. The behavior of Raa can be 
used to provide information on signatures of quark-gluon 
plasma formation [11 [H [2(| [H| ■ Pion parameterizations 
have been used for the proton-proton cross sections in 
these nuclear modification factors [H, [H HI . Given 
such widespread use of pion production cross section pa- 
rameterizations, we deem it worthwhile to test currently 
available parameterizations against new accelerator data. 
This is the aim of the present work. 

The outline of the paper is as follows. We summarize 
the paramaterizations studied previously @ making note 
of necessary variable transformations needed to describe 
theNA49 data set. We then compare these param- 
eterizations to the new data for pp reactions. Therefore 
the present work is a continuation of the previous paper 
[9j but applied to the new data. The parameterizations 
which are able to give a reasonable fit to the pp data are 
then scaled and compared to the pC data. 



II. DIFFERENTIAL CROSS SECTIONS 



A. Review of Kinematics 



Consider the inclusive reaction 
a + b — > c + X 



(1) 
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where c is the produced particle of interest and X is any- 
thing. Throughout this paper we assume that all vari- 
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ables, such as all momenta, are evaluated in the center 
of momentum (CM) frame. The momentum of particle 
c is denoted as p, and supposing that it comes out at 
angle 9 to the beam direction, then the longitudinal and 
transverse components of momentum are 



p z = p cos 9 
Pt = psm.6 

Feynman used a scaled variable instead of p 7 itself 
23l 2 4 25 11 . The Feynman scaling variable is [24L I26L 

H Ha Ha IP 



(2) 
(3) 



Xp 



Pz 



Pz 



(4) 



where p z is the longitudinal momentum of the produced 
meson in the CM frame, and p z max is the maximum mo- 
mentum of the produced meson given by (27l . |28| . [3lj 



\(s,m c ,m x ) 



4.s 



with 



A(s, m,, nij) = (s — mf — m j) 2 ~ 4m|m| 
Note that 



Pz 



Pn 



(5) 



(6) 



(7) 



Nagamiya and Gyulassy [28fl point out that if c is a boson 
with zero baryon number, then 



m x = mA + ms 



(8) 



in agreement with the p z max formulas of Nagamiya and 
Gyulassy [28| and Cassing [3l|]. The Fey nman scaling 
variable approaches the limiting value [29j 



X F 



2p z 



as s — > oo 



(9) 



Also it is obviously bounded in the following manner [24[ 

- 1< x F < 1 (10) 

Sets of variables that are often used are either (p, 9) or 
(Pz,Pt) Writing 



Xp 



\(s,m c ,m x ) 



4.s 



(11) 



shows that another useful and common variable set is 
(xf,Pt), which is used by the NA49 collaboration P, 0| 
in presenting their data. These variables are also used 
throughout the present work. Rapidity is defined as 



1, (E + Pz 

2 \E-p z 



(12) 



so that 



E = mx cosh y 
p z = mx sinh y 

where the transverse mass is defined through 



Pt 



E 1 



Pz 



(13) 
(14) 



(15) 



with m as the mass of the produced particle c.This gives 
yet another useful variable set (y,pr)- In the work below 
it will be necessary to write the rapidity in terms of the 
Feynman scaling variable as 



1, (y/x 
V = 77 log 1 



m T/Pz max + X F 



Xf 



For massless particles, E — p, so that y becomes 
cos 9 



7} 



1, 1 

— log 

2 6 1 - cos 9 



log tan ; 



(16) 



(17) 



This is called the pseudorapidity and is a good approxi- 
mation to the rapidity for particles moving near the speed 
of light. Because the pseudorapdity depends only on an- 
gle it can be used as an angular variable. Wong (2?| 
provides useful formulas involving rj and also gives ex- 
pressions relating y to r\ for slower-than-light particles. 



B. Parameterizations 

Blattnig et al Q did a study of the various parameteri- 
zations available for inclusive pion production in proton- 
proton collisions. They concluded that the Badhwar pa- 
rameterization (To| worked the best for charged pion pro- 
duction. However other parameterizations [jl, [H, S, 
135] will be reviewed again to see which works best for the 
new experimental data. The NA49 data set [l], Q uses 
the variables (xf,Pt), whereas some of the parameteri- 
zations below are written in terms of other variables sets. 
These will need to be converted to (xf,Pt)- 



1. Badhwar parameterization 

This parameterization [lol ] gives the Lorentz-invariant 
differential cross section as 



E 



d 3 a 



A 



d 3 p (1 + Aml/s) r 



(1 - x) q exp[ 



-Bp T 



1 + Ami/ s 



(18) 



where m p is the proton mass, y/s is the total energy 
in the center of momentum (CM) frame, and px is the 
transverse momentum of the produced meson in the CM 
frame. The other terms are given by 



1/2 



(19) 



3 



where it is assumed that the variables appearing in xf 

!S 

Ci + C 2 pt + c 3P 2 T 



Ellis parameterization 



are in the CM frame. Badhwar writes xt = xp. Also 



(20) 



The constants are listed in Table 1. The Badhwar 
variables are (xf , Pt) .which are the variables used in 
the NA49 data set Q, Q , so that no variable conversion 
is necessary. 

TABLE 1. Constants for the Badhwar parameterization. 
Units for A, C 2 and C 3 are mb/GeV 2 , GeV" 1 and GeV" 2 
respectively, and other constants are dimensionless. 



Particle 


.4 


B 


r 




c 2 


c 3 


7T+ 


153 


5.55 


1 


5.3667 


-3.5 


0.8334 


7T~ 


127 


5.3 


3 


7.0334 


-4.5 


1.667 



The Ellis [33j parameterization is 
d 3 p 



= a( p 2 t + m z )-»i*(i - XT y 



(23) 



where A is an overall normalization fitted to be A = 13 
in reference @ and xt = Pr/Pmax ~ 2pr/y/s- The 
same value of A is used in the present work. The 
other constants are listed in Table 3. The Ellis 
parameterization is independent of the emission angle 
9, and so does not carry any dependence onp z , xf-, y etc. 

TABLE 3. Constants for the Ellis parameterization. 



Particle 


N 


M 2 


F 


7T+ 


7.70 


0.74 


11.0 




7.78 


0.79 


11.9 



2. Alper parameterization 

The Alper [32| parameterization used in reference @ 
was 

E< %>^ A ex P(~ B P T + °Pt) e M-Dy 2 ) (21) 

where A, B, C, and D are constants that depend on the 
value of y/s. A more general formula is [32| 



E 



fa 

d 3 p 



A\ exp(-BpT) exp(— Dy 2 ) 

(1 -PT/Pbcam)" 1 



"An 



{p 2 T + M 2 ) n 



(22) 



The constants are listed in Table 2. The Alper variables 
are (y^pr)- To change to the variables (xf,Pt), we 
convert the rapidity in (|2"2"]) to if using (JTJ 



TABLE 2. Constants for the Alper parameterization. 



Particle A x B 



D Ay 



M 



4- Mokhov parameterization 



The Mokhov [3J] parameterization is 

B 



where 



P 



Pu 



exp 



P 



C^fs 



Vi{p T )V 2 {pj 



Vxipr) = (l-D)cxp(-Ep 2 T ) + Dexp(~Fp 2 T ) 

forpT < 0.933 GeV 

0.2625 



{p 2 T + 0.87) 4 



for p T > 0.933 GeV 



and 



V 2 (pt) 



0.7363 exp(0.875p r ) forp T < 0.35 GeV 
1 for pt > 0.35 GeV 



The constants are listed in Table 4. Using p = \/p 2 + p T , 
gives the Mokhov variables (p z ,Pt) which are trans- 
formed to {xf,Pt) using (fTTj) . 

TABLE 4. Constants for the Mokhov parameterization. 



Particle 



.4 



B 



C 



D 



E 



F 



tt+ 210 7.58 0.20 10.7 1.03 10.9 4.0 n + 60 1 1 9 18 3 12 2 7 

7T- 205 7.44 0.21 12.8 1.08 13.1 4.0 n - 5L2 2.6 0A7 0.3 12 2.7 



4 




FIG. 1: 7r + production in pp collisions. Data from reference 
[lj is plotted against the parameterization of Badhwar [Iol |. 
The values of xf from top to bottom are 0.0, 0.01, 0.02, 0.025, 
0.03, 0.04, 0.05, 0.06, 0.07, 0.075, 0.08, 0.1, 0.12, 0.15, 0.2, 
0.25, 0.3, 0.35, 0.45, 0.55. Following QJ, data and lines are 
multiplied successively by 0.5 to allow for a better separation. 



5. Carey parameterization 

The Carey [35j parameterization, which only applies 
to 7r~ , is 

E % { ^ ] = N{p2r + °- 86 )~ 4 ' 5 ( 1 - ( 24 ) 

where N is an overall normalization fitted to be N = 13 
in reference 0] and xr = p/p max ~ 2p/y/s. The same 
value of N is used in the present work. The Carey vari- 
ables are (p z ,Pt)- To change to the variables (xf,Pt), 

we USe XR = y/ X 2 p + p T / 



III. COMPARISON TO DATA 

The above parameterizations are compared to the new 
experimental data in Figs. 1-13. 



FIG. 2: ir + production in pp collisions. Data from reference [H 
is plotted against the parameterization of Alper [32|. Values 
of xf and data and line multiplication are the same as Fig.l. 



A. Comparison to pp data 

Positive pion production in pp reactions is shown in 
Figs. 1-4. In agreement with the conclusions of the 
NA49 collaboration [l[ , none of these arithmetic param- 
eterizations is able to account for all the fine structure 
seen in the data. The Badhwar, Alper and Ellis parame- 
terizations are unable to reproduce the shape of the dif- 
ferential cross section at low pr- The Alper and Ellis 
parameterizations are unable to account for the data at 
larger values of xf- However, on the positive side, the 
Mokhov parameterization (Fig. 4) does give a reasonable 
description of the shape at low px and the Badhwar pa- 
rameterization (Fig. 1) gives a reasonable fit to the data 
over &\Ixf values for px > 0.3 GeV. In general, the Alper 
and Ellis fits are not good, but the Badhwar and Mokhov 
fits are moderately good, but certainly not precise. 

Negative pion production in pp reactions is shown in 
Figs. 5-9. Again, in agreement with the conclusion of 
reference [l[ , none of these arithmetic parameterizations 
is able to account for all the fine structure seen in the 
data. The conclusions are the same as for positive pions 
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FIG. 3: 7r production in pp collisions. Data from reference 
[lj is plotted against the parameterization of Ellis [U . Values 
of xf and data and line multiplication are the same as Fig.l. 



FIG. 4: 7r production in pp collisions. Data from reference 
[lj is plotted against the parameterization of Mokhov [34| . 
Values of if and data and line multiplication are the same 
as Fig. 1. 



production, except that the Badhwar parameterization 
now fails for larger values of xp. The Carey parameteri- 
zation is reasonable except for small values of xp at low 
Pt- In general, the Alper and Ellis fits are not good, 
but the Badhwar, Mokhov and Carey fits are moderately 
good, but not precise. 



B. Comparison to pC data 

Although the primary emphasis of this paper is on the 
pp reactions, nevertheless it is of interest to see if the pp 
parameterizations are able to scale to fit the new pC data 
[2J, although a precise fit is not expected. We will only 
consider the parameterizations of Badhwar and Mokhov 
because they provided the best fits to the pp data. 

One should note that in comparing to the pC data 
we will assume the same cross section for proton-proton 
and proton-neutron scattering. However, as shown in 
the paper by Pawlowski and Szczurek [36| , this is known 
not to be true in general for the energy region of the 
new NA49 data. We justify our use of the same cross 



sections by noting the following. The data of Pawlowski 
and Szczurek [36| (see their Fig. 2), show that these 
cross sections are significantly different only for xp < 
—0.05. For xf > —0.05 the cross sections are very similar 
and in fact for xp > 0, Pawlowski and Szczurek (36j 

themselves state that a^ p « a^ n . In our comparisons, 
we show 25 different values of xp. Only two of these, 
namely xp ~ —0.01 and xp = —0.075 are such that 
there is a significant difference between the proton-proton 
and proton-neutron cross section. For 7r + production the 
proton-proton cross section is about twice that of the 
proton- neutron cross section and vice- versa for tt~ [36| . 
This means that only for xp — —0.01 and xp — —0.075 
our comparisons to data will be slightly worse for tt + and 
slightly better for 7r~. We now discuss various models. 

The wounded nucleon model was introduced by Bialas, 
Bleszynski and Czyz [37| . The number of wounded nu- 
cleons is simply the number of participants involved in 
the reaction. The main assumption is that the parti- 
cle multiplicty is proportional to the number of wounded 
nucleons. Consider pA reactions. The incident proton in- 
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FIG. 5: ir~ production in pp collisions. Data from reference [l| 
is plotted against the parameterization of Badhwar [hj. The 
values of if from top to bottom are 0.0, 0.01, 0.02, 0.025, 
0.03, 0.04, 0.05, 0.06, 0.07, 0.075, 0.08, 0.1, 0.12, 0.15, 0.2, 
0.25, 0.3, 0.35, 0.45, 0.55, 0.65, 0.75, 0.85. Following [fl], data 
and lines are multiplied successively by 0.5 up to if = 0.35 
and by 0.75 for xf > 0.45 to allow for a better separation. 



teracts with v nucleons in the target nucleus [38|, which 
is some fraction of the total nucleon number A, and 
is determined from collision geometry and the hadron- 
nucleon cross section (38|. In a pA reaction, the number 
of participants is (1 + v). Let N p a be the multiplic- 
ity of particles of interest produced in the pA reaction, 
and let N pp be the particle multiplicity in the pp reac- 
tion. For example this could be the number of pions 
produced in a reaction. The multiplicities are related by 
N P A = 4(1 + v)N pp The | factor is there because the 
"nucleon-nucleon interaction requires two wounded nu- 
cleons" [39| . The following comes from reference (39| . 
For a nucleus-nucleus collision let A be the number of 
nucleons in the projectile and let B be the number of 
nucleons in the target. Then [3!| Nab — \WabN pp 
where Wab is the number of participant nucleons, i.e. 
the "number of nucleons that have interacted at least 
once" H. This is given by [H W AB = A^ + B^ 
where <j p a and a P B are the proton-nucleus inelastic cross 



FIG. 6: 7r~ production in pp collisions. Data from reference [l| 
is plotted against the parameterization of Alper [3^|. Values 
of if and data and line multiplication are the same as Fig. 5. 



sections and gab is the nucleus-nucleus cross section 
given in reference [39]. If the projectile nucleus is ac- 
tually a proton then W vB = ^-+B^ = l + B^- or 

l + v= W vA = I + AZvl showing that \M v = A^-. In 
pA or AA collisions then the particle multiplicity scales 
with the number of wounded nucleons, which is often 
calculated with Glauber theory (4lj |. 

Instead of wounded nucleon or participant scaling, one 
can have binary scaling. The total number of produced 
particles in nuclear collisions comes from both soft and 
hard processes. Soft processes should scale with the num- 
ber of participants and hard processes should scale [42j 
with the average number of binary nucleon-nucleon colli- 
sions denoted as iV co n. Reference [43| contains a nice in- 
troduction to hard processes. For example, the PHENIX 
experiment [431 ] focuses on detecting large transverse mo- 
mentum pt particles that arise from the early stages 
of relativistic heavy ion collisions. In the early stage 
nucleon-nucleon collisions cause jet production resulting 
from hard parton collisions. The jets subsequently decay 
into high momentum hadrons with p? ~ 2 Gev [431 ] . The 
quark-gluon plasma (QGP) forms at a later time in the 
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FIG. 7: 7r~ production in pp collisions. Data from reference 
[H is plotted against the parameterization of Ellis [33|] ■ Values 
of xf and data and line multiplication are the same as Fig. 5. 



FIG. 8: 7r~ production in pp collisions. Data from reference 
[lj is plotted against the parameterization of Mohov [34|] . Val- 
ues of xf and data and line multiplication are the same as 
Fig.5. 



collision and the early scattered partons move through 
the QGP region leading to the phenomenon known as jet 
quenching which is signified by a "depletion in the yield of 
high pt hadrons" [43(] . Observation of this jet quenching 
is therefore a "potential signature for QGP formation." 

The nuclear modification factor is a measure of nuclear 
effects and has been discussed by many authors [H, 52, 
IH, 0, HE El] ■ For proton- nucleus reactions, 



R-oA = 



Ed 3 a p A/d 3 p 



TVcoii Ed 3 a pp /d 3 p 
For nucleus- nucleus reactions, 



Rab = 



Ed 3 a A B/d 3 p 
N co \\ Ed 3 a N N/d 3 p 



(25) 



(26) 



where NN refers to the nucleon-nucleon cross section. 
Instead of d 3 p, other variables can be used such as xf, Pt, 
rapidity, psuedorapidity, etc. The nuclear modification 
factor R = 1 in the absence of nuclear effects, i.e. if the 
nuclear cross section is just an incoherent superposition 
of nucleon-nucleon collisions. For low px < 2 GeV it has 
been found I43| that R < 1 and this is due to the fact 



that the reactions scales with the number of participants 
(participant scaling), rather than the number of binary 
collisions (43|. For pr > 2 GeV "particle production in 
pA collisions is enhanced compared to binary scaling" 
[43j |. which is the Cronin effect. 

Again consider participant scaling versus binary scal- 
ing. The nucleus-nucleus cross section could scale either 
as a function of the number of participants (often called 
wounded nucleons) or as a function of the number of bi- 
nary nucleon-nucleon collisions. Binary scaling would in- 
dicate no collective effects whatsoever. However if there 
are some collective effects then these should manifest 
themselves by scaling with the number of participants. 
Obviously binary collisions are harder and lead to high 
Pt processes such as jet production (from individual par- 
ton collisions) and heavy flavor production. Participant 
effects are generally softer processes at smaller pr such as 
soft hadron production, transverse energy flow and other 
collective phenomena. Both the number of participants 
Apart and the number of binary collisions N co \\ is large for 
small impact parameter (central collisions) and drops off 
smoothly at larger impact parameters (peripheral colli- 
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FIG. 9: 7r~ production in pp collisions. Data from reference 
[lj is plotted against the parameterization of Carey [35j]. Val- 
ues of if and data and line multiplication are the same as 
Fig.5. 

sions). Also we always have iV co n > iV par t, because there 
can be many re-scatterings among the participants. 

For a long time it has been known that energetic 
particle production in proton-nucleus collisions increases 
faster than the number of binary nucleon-nucleon colli- 
sions p6| . In other words, particle production in nuclear 
collisions is enhanced compared to binary scaling [431 ] . 
This is known as the Cronin effect. Various physical 
mechanisms can contribute to the Cronin effect such as 
multiple parton scattering in the initial stage of the col- 
lision [42j| . Fermi motion |44l| . etc. Reference [46[ points 
out that the "the cause of the Cronin effect and its species 
dependence are not yet completely understood", where 
species dependance refers to the fact that the size of the 
effect varies for different particles produced. A contribu- 
tion from final state interactions is also possible [46[ . Ref- 
erence [46] provides a very nice summary of the Cronin 
effect. The pA cross section is parameterized as (-461. |47| 

E %^ A ) = E %(^) Aa[PT) ( 27 ) 
where a > 1 indicates nuclear collective effects. Adler et 



FIG. 10: 7r + production in pC collisions. Data from reference 
[2] is plotted against the parameterization of Badhwar pTj ] 
multiplied by a best fit factor of 12 0,9 . The values of from 
top to bottom are -0.1, -0.075, -0.05, -0.04, -0.03, -0.025, -0.02, 
-0.01, 0.0, 0.01, 0.02, 0.025, 0.03, 0.04, 0.05, 0.06, 0.075, 0.1, 
0.125, 0.15, 0.2, 0.25, 0.3, 0.4, 0.5. Following data and 
lines are multiplied successively by 0.5 to allow for a better 
separation. 



al [46( state that "the enhancement depends on the mo- 
mentum and the type of particle produced, with protons 
and antiprotons exhibiting a much larger enhancment 
than pions and kaons at px > 2—3 GeV." They also point 
out that at y/s = 27 .4 GeV, the enhancement reaches a 
maximum at pt = 4.5 GeV with ~ a+ = 1.1 and 
that at the same momentum a p = 1.3 for protons [46l |. 

In Figs. 10 - 13 we have simply multiplied the pp pa- 
rameterizations of Badhwar and Mokhov by a constant 
Cronin enhancement factor, 12°, and compared them to 
the pC data from NA49 0. In our previous pp fits, the 
Badhwar parameterization worked best for high px > 0.5 
GeV and the Mokhov parameterization worked best for 
low pt < 0.5 GeV. Therefore we do not expect the Bad- 
hwar parameterization to fit low pr data for the pC reac- 
tion and we do not expect the Mokhov parameterization 
to fit high pt data for the pC reaction. We varied the 
value of a to find the best possible fit to the pC data. 



100 




100 



(Badhwar ) 





FIG. 11: 7r production in pC collisions. Data from reference 
0| is plotted against the parameterization of Mokhov [34| 
multplied by a best fit factor of 12 0,9 . Values of if and data 
and line multiplication are the same as Fig. 10. 



FIG. 12: 7r~ production in pC collisions. Data from reference 
is plotted against the parameterization of Badhwar fiol ] 
multplied by a best fit factor of 12 0,9 . Values of if and data 
and line multiplication are the same as Fig. 10. 



Figs. 10-11 show the Badhwar fit to the high px data 
with a best value of a = 0.9 and Figs. 12 - 13 show the 
Mokhov fit to the low px data also with an independent 
best value of a — 0.9. 



IV. CONCLUSIONS 

The NA49 collaboration [J, has provided new pre- 
cise data on pion production in pp and pC reactions at 
a beam momentum of 158 GeV. Although a precise fit 
is not expected, nevertheless it is of interest to compare 
currently available arithmetic parameterizations to the 
new data. Let us emphasize that we are not suggesting 
that arithmetic parameterizations are able to give a com- 
plete account of the new data et al [H, 0] . The numerical 
interpolation developed by the NA49 collaboration [J, Q 
is far superior. The aim of this paper has rather been to 
see how well some parameterizations describe the data. 

We conclude that the Alper and Ellis parameteriza- 
tions do not fit the pp data very well, and should not be 



used in this energy region. The Carey parameterization 
for 7r~ works better but underpredicts the data at low pt 
for small values of Xp and the predictions at high px are 
only moderately good. The Badhwar parameterization 
for tt^~ works well for high values of px, but does a poor 
job at low pt, whereas the Mokhov parameterization is 
the other way around. It works best for low pr but not 
well for high px- Note also that the Badhwar parameter- 
ization for 7r~ docs not work well for large values of xp, 
whereas the Mokhov parameterization works fine in this 
region. Regarding the pp parameterizations, we conclude 
that for low px < 0.5 GeV, it is best to use the Mokhov 
parameterization. For high px > 0.5 GeV, is best to use 
the Badhwar parameterization, except for large values of 
xp > 0.45, where it is better to use Mokhov. 

The Badhwar and Mokhov parameterizations gave 
good fits in certain px ranges and we scaled them to 
the pC data. In the Cronin effect, the scaling factor is 
A a (PT) _ -yy e f 0U nd that the Badhwar parameterization 
gave the best fit for a w 0.9, and the Mokhov parame- 
terization also gave the best fit for a sa 0.9. As discussed 
above, the Badhwar parameterization works best for high 



10 



100 



p + C -> 7T + X (Mokhov) 





U 




cn 


qui 


> 




0) 




U 







a 




ro 




X! 


H 



0.01 



U 0.0001 



1. x 10" 



1. x 10" 8 




Pt and the Mokhov parameterization works best for low 
Pt- Therefore we conclude that the Cronin enhancement 
factor for the pC reactions is a(pr) ~ 0.9 in both the 
high px > 0.5 GeV region and low px < 0.5 GeV region. 
These results indicate the absence of nuclear collective 
effects [4(| for this pC reaction Q ■ 
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FIG. 13: 7T~ production in pC collisions. Data from reference 
is plotted against the parameterization of Mokhov [34j] 
multiplied by a best fit factor of 12 ' 9 . Values of if and data 
and line multiplication are the same as Fig. 10. 
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